Velocity-force characteristics of an interface driven through a periodic potential 
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We study creep dynamics of a two-dimensional interface driven through a periodic potential using 
dynamic renormalization group methods. We find that the nature of weak-drive transport depends 
qualitatively on whether the temperature T is above or below the equifibrium roughening transition 
temperature Tc. Above Tc, the velocity-force characteristics is Ohmic, with linear mobility exhibiting 
a jump discontinuity across the transition. For T < Tc, the transport is highly nonlinear, exhibiting 
an interesting crossover in temperature and weak external force F. For intermediate drive, F > F^, 
we find near a power-law velocity-force characteristics v{F) ~ F"^ , with a — 1 oc t, and well- 
below Tc, v(F) ~ e-iP*/Ff\ with i = {I - T/Tc). In the limit of vanishing drive {F < F,) the 



velocity-force characteristics crosses over to v{F) 
I. INTRODUCTION 

The problem of elastic media pinned by an exter- 
nal potential provides a unifying framework for un- 
derstanding a large number of condensed matjtj^r phe- 



-(Fo/F) 



, and is controlled by soliton nucleation. 



nomena, such as for example surface-growtlu'D, non- 
linear transport|-.ia anisotropic metelaj, dissipation in 
supercaoductorsoQ, Wigner crystalsu, earthquakeaa and 
frictionllJ. Randomly-pinned elastic systems are also toy 
models_for considerably more complicated problems of 
glassesEJ. Much of the recent interest in -such prob- 
lems has been rekindled by the discoverytJ of high- 
temperature superconductors (HTSC), and efforts to un- 
derstand the nature of their B — T phase diagram and 
dissipation controlled Jaw, siatics and dynamics of elastic 
arrays of vortex linesE^cm. A combination of thermal 
fluctuations, pinning, and external drive leads to a wide 
range of new and interesting collective phenomena that 
is common to many physical realizations of elastic media. 

There has been considerable progress in undep 
ing the static properties of pinned elastic media 
Much of the recent interest has therefore shifted to dy- 
namics, with current focus aiiiaonequilibrium, driven dy- 
namics of these rich systemsLTLl. Once the elastic medium 
is driven, however, many new questions arise, such as 
the governing nonequilibrium equation of motion, phase 
classification and stability, nature of the corresponding 
phase transi tions , and the resulting nonequilibrium phase 
diagramfl"l3'EIrEj. Among these many challenging ques- 
tions, the velocity {v) - force {F) characteristics (the 
IV curve, in the context of superconductors and charge 
density waves (CDW's)) is the observable that is most 
directly accessible experimentallyllj and is therefore of 
considerable theoretical interest. 

Despite ofj^onsiderable richness of many aspects of the 
driven stateOi3't3~E3, at large drives the velocity-force 
characteristics of a uniformly sliding medium approaches 
Ohmic form with deviations 5v that can be computed 
perturbativelyEil in the ratio 5v/v. At zero temperature, 
if elasticE3, the medium is pinned for drives smaller than 



a critical T = value Fc , and undergoes a nonequilibrium 
depinning transition to a sliding state, with v ~ \F — Fc\^ 
playing the role analogous to an order parameter dis- 
tinguishing pinned and sliding phasesci Finite temper- 
ature rounds the depinning transitioro, allowing acti- 
vated creep motion of the elastic solid even for drive far 
below Fc- 

Much of the insight about this highly nontrivial creep 
regime that is the focus of our work come.s .from-.a scal- 
ing theory of depinned droplet nucleationE£rE3c3. This 
approach generically predicts collectively pinned elastic 
media to exhibit a highly nonlinear v{F), with a van- 
ishing linear mobility, corresponding to transport acti- 
vated over barriers that diverge with system size and 
vanishing drive. Recently, in the case of random pin- 
ning, these scaling predictions have been put on firmer 
ground through a detailed dynamic functional cepormal- 
ization group (DFRG) calculations of i;(F)l23'Ej, which 
indeed predict v{F) ~ e~^/^^, with a universal fi expo- 
nent. However, in the case of random pinning a number 
of technical problems with DFRG remain, precluding a 
fully controlled analysiSc3. 

It turns out, however|^hat many problems of interest, 
such as surface-gmwtho, 2D colloidal crystals in peri- 
odic pptentialsE3~L3, and vortices pinned by artificial dot 
arraygSa or by intrinsic pinning in e.g., HTSC, involve 
the considerably simpler but still nontrivial problem of 
periodic pinning. In addition to addressing numerous in- 
teresting physical problems, study of motion in a periodic 
potential provides a nice laboratory to explore new cal- 
culational methods. 

A driven sine- Gordon model is the simplest descrip- 
tion of such periodic pinning problem, with an impor- 
tant simplifying feature of absence of topological defects 
such as dislocations or phase slips that can be impor- 
tant for undetatanding the dynamics of vortex arrays 
and CDW''a3'c2l. Directly applicable to crystal growth 
phenomena, tius-jpodel has been extensively studied in 
the literatureEirE2l. In equilibrium, among many other 
things, it describes the famous crystal surface roughen- 
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ing transition from the low-temperature smooth phase 
with bounded surface roughness to the high-temperature 
rough phase with logarithmic height correlations. 

One of the many interesting questions that naturally 
arises is: What are the signatures of the roughening tran- 
sition in the driven transport ? More specifically here we 
are interested in qualitative differences (if any) in v{F) 
above and below Tc- The equilibrium limit of this ques- 
tion was addressed in the classic analysis of the equi- 
librium dy|iiamics by Chui and Weekscj and by Nozieres 
and Gallctca. By extending the standard analysis to equi- 
librium dynamics they found the vanishing of the linear 
mobility below the roughening transition, consistent with 
the static picture of the smooth phase where the interface 
is pinned by the periodic potential. 

In the presence of an external drive, the sine-Gordon 
model was studied sometime ago b5jJIwa, Kardar, and 
Paczuskitil and by Rost and Spohnca. Although these 
works led to considerable progress, computing the rcnor- 
malization group (RG) flow equations for model parame- 
ters at nonzero F, they concentratedjnainly on the influ- 
ence of Kardar-Parisi-Zhang (KPZ)E3 nonlinearities (not 
considered in previous studiesH) important at strong ex- 
ternal drive, but said little about the actual v{F) char- 
acteristics of the driven interface in the _F — > creep 
regime. The driven sine- Gordon model has also been con- 
sidered by Blatter et ala in a complementary approach 
via a high velocity perturbative expansion for v{F). Con- 
sistent with Refs. ^ and these last authors found 
that while the correction {6v/v) remained flnite above 
Tc, it diverged below Tc, thereby suggesting nontrivial 
transport changes across the roughening transition, but 
leaving the form of v{F) in the creep regime an open 
problem. Of course, because at flnite temperature the 
interface can move for any flnite drive F, at sufficiently 
long scales the periodic potential is averaged away at 
both low and high temperatures, thereby leading to the 
rounding of the roughening transition itself. Neverthe- 
less, we expect that the velocity-force characteristics in 
the creep regime is controlled by the equilibrium physics 
and precise qualitative distinction of v{F) in the rough 
and smooth phases should existEJ. 

In this paper, our goal is to understand in detail the 
physical consequences of the divergences found in the 
high- velocity perturbative expansion and in particular to 
compute the creep velocity-force characteristics in both 
phases and,acrpss the roughening transition, utilizing dy- 
namic RGO'ErO. Consistent with perturbative analysis, 
we find that the nature of transport depends qualitatively 
on whether the temperature is above or below the equi- 
librium roughening transition temperature Tc- Above 
Tc, the velocity-force characteristics is Ohmic, with the 
mobility remaining finite for T — > T^. In contrast, for 
T < Tc, we find that the linear mobility vanishes on long 
length scales, and therefore exhibits a nonuni versa 1. jump 
discontinuity across the roughening transitiorOO'u. In 
the smooth phase, the transport is a strongly nonlin- 
ear function of applied force, showing a rich universal 



crossover in temperature and applied force (Fig. |^). 
At an intermediate drive F > F^,{g,T), larger than the 
pinning (g)- and temperature-dependent strong-coupling 
crossover force 



F*i~g,T) 



T —*T- 
T<^Tc 



(1.1) 



the velocity-force characteristics strongly depends on the 
level of proximity to Tc, with: 



v{F) 



T^T- 
T<^Tc 



(1.2) 



where t = (I — T/Tc), and bi and &2 are nonuniver- 
sal constants of order unity. For sufflciently low drive 
F < F^{g, T), the motion is instead always via activated 
soliton creep, with the velocity-force characteristics cross- 
ing over to 



v{F) 



-Fo/F 



F<F„ 



(1.3) 



with To another characteristic force that will be deflned 



below, Eq.(^.17D. 

This paper is organized as follows. We introduce the 
driven sine-Gordon model in Section || and analyze it in 
Section III using simple perturbation theory in the 




FIG. 1. Typical velocity-force characteristics of a driven 
interface in a periodic potential. At zero temperature, the 
interface remains pinned {v = 0) until F reaches the criti- 
cal force Fc = pg. At finite temperatures, we find that the 
near-equilibrium response of the interface to a small (T <^ Fc) 
driving force depends on whether T is above or below the 
roughening temperature Tc- For T > Tc, the velocity- force 
characteristics is Ohmic {v{F) ~ F) down to T = 0, while 
for T < Tc and forces smaller than a characteristic force T*, 
the characteristics is strongly non-linear, v ~ exp ( — Fq/F^, 
creep motion via activation over barriers that diverge in a 
vanishing drive limit. 
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pinning potential strength. While for weak pinning this 
computation is convergent for T > Tc, it fails for ar- 
bitrarily weak pinning in the smooth phase. In section 
IV we employ dynamic RG techniques to make sense of 
these divergences, and in Section ^ we use these results 
to compute v{F) through the roughening transition. We 
conclude in Section VI with a summary of the results and 
a discussion of open problems and future directions. 



II. DRIVEN SINE-GORDON MODEL 

In equilibrium a two-dimensional sine-Gordon model 
of an elastic interface is described by a Hamiltonian 



H ^ dr 



g cos (p/i(r)) 



(2.1) 



where r is a two dimensional vector in the (xy) plane, 
h{r) is the height of the interface above the {xy) plane 
(taken to be along the z direction in the embedding 
space) at location r, K is the interfacial surface tension, 
g is the pinning strength and d = 271 /p is the period of 
the potential. In the context of a crystalline surface, with 
H characterizing its equilibrium roughness, the periodic 
pinning potential softly encodes lattice periodicity of the 
bulk crystal, corresponding to the ft, — > h + d a symme- 
try of the surface energy, with d being the crystal lattice 
constant perpendicular to the interface. 

In the absence of any additional conservation laws, long 
scale equilibrium dynamics can be described by a simple, 
relaxational (model A) Langevin equation 



jdth 



5H 



+ C(r,t) , 



Sh{r,t) 

KV^h{r, t) - pg sin {ph{r, t)) + C(r, t) , (2.2) 



where 7 is the microscopic friction coefficient, and C(r, t) 
a zero-mean, Gaussian thermal noise describing the in- 
teraction of the system with the surrounding heat bath 
at temperature T, with 



{Cir,mr',t'))=2jTSir-r')S{t-t'), 



(2.3) 



in equilibrium imposed by the fluctuation-dissipation 
theorem (FDT), forbidding independent renormalization 
of T. 

The dynamic description of an interface driven by an 
external force F (in the context of crystal growth propor- 
tional to the difference between the chemical potentials 
of the solid and vapor phases) is substantially modified. 
In addition to the obvious addition of the driving force 
F on the right-hand side of Eq. ( |2.2| ) , nonequilibrium dy- 
namics permits the appearance of nonconservative forces 
(those not expressible as derivatives of H), the most iija^. 
portant of which is the famous Kardar-Parisi-Zhang£3 
(KPZ) (V/i)^ nonlinearity, allowed by the explicit break- 
ing by the drive of the z —z symmetry. An additional 



important effect of driving appears as the renormaliza- 
tion of "temperature" T, corresponding to the breakdown 
of the fluctuation-dissipation theorem, that is, the renor- 
malization of the friction coefficient 7 is independent of 
that of the variance of the noise ^. Even if these nonequi- 
librium effects are not recognized a p riori, they appear 
upon coarse-graining of eaiLation (2.2) as soon as the ex- 
ternal drive F is includecS'u. The pssulting nonequilib- 
rium equation of motion is given byEJ 

J dth^KV'^h+^{Vhf~pg sin {ph{r,t)) +F + C(r,t) . 

(2.4) 

Our goal here is to apply the machinery of the dynamic 
RG to compute the velocity-force v{F) characteristics for 
the above model, focusing on the nontrivial creep regime 
of the smooth phase, where naive perturbative expansion 
in the pinning potential g fails. 



III. DYNAMIC PERTURBATION THEORY 

It is instructive to first study the velocity-force charac- 
teristics through a simple perturbative expansion in the 
pinning potential g. Starting from equation (2.4), it is 
convenient to shift h{r,t) — vqI + u{r,t) with vq = F/7 
the unperturbed {g — X = 0) expression of the velocity. 
Averaging ( ^.4[ ) over thermal fluctuations, and ignoring 
the KPZ term, we find that the velocity v of the moving 
interface is given by 



V = (dth) (3.1a) 

^ ~ ~~ (si'^ (p"(i", t) + pvot)) , (3.1b) 

where we used the fact that (C(r,t)) = 0. We now let 

u{r,t) ^ uo{r,t) +Ug{r,t) , (3.2) 

where 



uo{r,t)^ J dr'dt' Rn{r-r',t-t')C{r\t') 



(3.3) 



is the thermal (nonintcracting) part of the interface dis- 
placement, 

Ug{r, t) =pg j dr'dt' i?o(r - r', i - t') 

{P^+pu,{v',t')) (3.4) 



X sm 



is the correction to u linear in the pinning potential 
strength g, and Ro{r - r',i - t') = 6{uQ{r,-l^/ 5F{y' ,t') 
is the response function of the free interfaceEj. Expand- 
ing equation ( 3.1b| ) in Ug, and avera ging oy er the thermal 
noise we find (see also Appendix VII B) 
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7 27 7 

vF 

X sin[— (t-t')]^o(r-r',t-i') , (3.5) 

7 

where Co(r - r',t - t') = ([uo(r,t) - uo(r',t')]^> is the 
connected correlation function of a free interface given 

by 

Co(r,i)c.^ ln[l + A2(r^ + ^)] . (3.6) 

Th e above velocity-force characteristics, equation 
( |3.5|) , is most easily evaluated at zero temperature where 
Co(r, i) = 0. In this limit, using 



i?o(r,t) 



e{t) 



(3.7) 



{9{t) is Heaviside's unit step function) inside Eq.(3.5) and 
integrating over the time variable t', we obtain 



F p^g^ 
7 27 



dr 



pF 



(3.8) 



In the above equations and througout the rest of this 
paper we use a shorthand notation for /c?q/(27r)^, 
and the superscript A — 211 /a is the ultra-violet cutoff 
set by the in-plane lattice constant a, generically distinct 
from the period d — I-k jp perpendicular to the interface. 
Performing the integration over the space variable r in 
the last equation, and using the resulting Dirac delta- 
function (27r)^(5(q) to complete the q integral, we find a 
T = 0, leading ord^p-^^^ pinning g) expression for the 
viF) characteristics 



7 V 2 ^ F'' 



(3.9) 



where F^ = pg is the zero-temperature critical force, in 
agreement with the condition Fc — m&yS^dV {h) j dh\ of 
disappearance of metastability iy(Ji) — —gcoa{ph) is the 
pinning potential). As is clear from this result for v{F), 
even at T = 0, the perturbative corrections are small 
for sufficiently large applied force F relative to the pin- 
ning force Fc (equivalently, for sufficiently weak pinning 
g at fixed F). In this fast moving regime, the metasta- 
bility is absent and pinning gives only a small correction 
to the motion with v{F) deviating only weakly from the 
pinning- free Ohmic response vq{F) = F/"f. It is reas- 
suring to note that, since at T = 0, only the q = 
mode contributes to the v{F), Eq. ( |3.9D agrg^ with the 
high-drive limit of the exact T — resultLjO for a sin- 
gle particle driven through a one-dimensional sinusoidal 
potential 



viF) 




F>Fc 



(3.10) 



This suggests that the v{F) characteristics of a driven 
interface should also exhibit a square-root cusp with an 
infinite slope at F = i^c- At T = the interface is strictly 
pinned for F < Fc. 

In contrast, at any finite temperature the interface 
moves for arbitrarily weak force and hence there is no 
sharp depinning transition. The perturbative expression 
for V, Eq. (3.5) can be readily evaluated by using the 
fiuctuation-dissipation relation 



e{t)dtCoir,t)=2TRoir,t) 



(3.11) 



obeyed by the equilibrium response and correlation func- 
tions. Usin g th is relation to eliminate Ro{r,t) from the 
rhs of Eq. (3.5) and integrating by parts over t' we find 



_ F r 

7 



^-^Jd^l cos{pFt/j)e- 



(3.12) 



Inserting into this last equation the expression of the cor- 
relation function Co(r, t) of a harmonic interface given in 
Eq. (Eel) leads to 



7 V 27r 



where we defined 



00 

dri dt 



cos{pFt/j) 



AttK' 



(3.13) 



(3.14) 



Taking the limit F — > in the above expression, and 
performing the time integration, we obtain 



^ui(v/F) = -\l 

F^O 7 



dr 



np^g^ 

KTA^ Jo (1 + A2r2)(2')-3)/2^ 

(3.15) 



We n ow observe that the integral on the rhs of equation 
( |3.15 ) behaves very differently depending on whether T 
is smaller or greater than 



Tco — 



8t:K 



p^ 



(3.16) 



For T > Tco, i.e., r] > 2, the integral in Eq. ( 3.15 ) is con- 
vergent, and leads to a finite (and for weak pinning g, to 
an arbitrarily small) correction to the linear friction coef- 
ficient 7(F = 0) = 1/ \imF^o{v/F). In strong contrast, 
for T < Tco [rj < 2) above integral diverges signalling the 
breakdown of the perturbation theory at small values of 
the external force F. 

Having established the breakdown of perturbation the- 
ory for T < Tco in the limit of vanishingly small forces, 
we now turn our attention to the full velocity-force char- 
acteristics at finite values of the external drive. Starting 
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from equation (3.13), and performing the integration over 
space variables, we obtain 



/ = 



pF 



(3.18) 



7 ^ 



1 



4 ? 

V 9 



V 



8i^2A4 - 1 



dT 



cos(2/r) 
(r + l)"-i 



(3.17) 



Performing the integralEl on the rhs of equation (3.17), 
we finally arrive at the following result for the effective 
friction coefficient j{f) of the driven interface (here 1F2 
is a generalized hypergeometric function) 



where the dimensionless force / is given by (henceforth, 
we shall use both F and / to designate the driving force 
on our interface) 



7(/)=7{l 



4 9 

P 9 



1 



SK'^A^ 77 - 1 
which has the following limiting behavior as / 



{2fr-'T{2 - r,) sin (2/ + -(^ - 1)) + _ 



Tj 3 — r] 



;-f)]y 



0, 



7(/^0)=7 



1 



4 2 

p 9 



8^2 A4 (2-?7)(r7-l) 



l-(2-77)r(2-77) (2/)''-2sin(-(7?-l 



(3.19) 



(3.20) 



As found above, inside the rough phase, T > T^q 
{rj > 2) and for sufficiently weak pinning, the perturba- 
tion theory remains valid at arbitrary /, simply display- 
ing crossover from a freely moving interface with "bare" 
mobility /ioo = I/7 at high drives to that with finitely 
suppressed low-drive mobility (as illustrated in figure H) : 



On the other hand, in agreement with Ref. 0, we find that 
in the "smooth" , low temperature T < Tc (77 < 2) phase, 
the behavior is stri kingl y different with the correction to 
vo{F) = F/7, Eq. (3.20), diverging and the perturbative 
approach failing as / is reduced below a characteristic 
force 



M9,T) 



S/f^A" (r)-l) 



sin(f(r;-l)) 



1 



^ (2-r,)(r,-l) 



1/(2-^) 



(3.22) 
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FIG. 2. Effective friction coefficient y{f) of the driven in- 
terface for rj = 2.2 and {p'^g^ /8K'^ A*) = 0.02. As / ^ 0, 7(/) 
remains finite, in agreement with equation (3.21). 
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FIG. 3. Effective friction 7(7) for rj = 1.8 and 
(p^ffVS-K^^A*) = 0.1. 7(/) diverges at /, ^ 0.022 (dashed 
line), indicating the failure of perturbation theory at small 
drives < / < /*. 



Asr->r-o (r;^2-), 

MT^To) ~ i exp( 



(3.23) 



showing that the regime of forces < / < /* where 
perturbation theory fails becomes exponentially small as 
Tco is approached from below. The unbounded growth of 
the perturbative friction coefficient as the external drive 
/ approaches /* from above (see Fig. ^) suggests that 
the interface in the low temperature, smoatk. phase is 
characterized by a vanishing linear mobility.Cj'Q 

Although we will study this in more detail in follow- 
ing sections, already at this stage we can see a physical 
interpretation of this divergence. Perturbation theory in 
the pinning potential fails because even for an arbitrarily 
weak pinning g, on sufhciently long scales greater than 
^ (computed in Section iV), the periodic potential (for 



small h acting like a "mass", hgp^h?) necessarily dom- 



inates over the elastic energy density (V/i) . Since 
(as is quite clear from the equation of motion, Eq. (^^)) 
the applied force F dominates the elastic force on scales 
longer than 
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6 



[pF J ■ 



(3.24) 



a sufficiently weak; force, F < F^,, probes the interface 
on length scales longer than ^ and thereby leads to the 
breakdown of perturbation theory about the harmonic 
interface. Hence, although quite instructive, the per- 
turbation theory fails to make predictions for v{F) or 
any other dynamic quantity in the smooth phase at suffi- 
ciently low drive F < F,, and a nonperturbative approach 
is necessary. 



IV. DYNAMIC RENORMALIZATION GROUP 

Armed with the above discussion, we are now well- 
equipped to use dynamic RG analysis to make physical 
sense of these perturbative divergences, with the main 
goal being the calculation of v{F) in the smooth T < Tc 
phase for weak drive F < F^. It is convenient to perform 
this analysis in the frame co-moving with the "bare" ve- 
locity vq = F/"f corresponding to the change of the dy- 
namic fields to u(r, t) = h{r, t) — Ft/^, which obeys 

A 2 pF 
7 dtu = KV'^u + - iyu) + pg sin(pu + ^t) + C{y, t) . 
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(4.1) 



Taking the nonlinear terms in the above equation as a 
small perturbation, the equation of motion, can be di- 
rectly expanded in these nonlinearitiea^J^LJ leading to 
renormalization group recursion relations for model pa- 
rameters. An equivalent but more convenient formula- 
tion ia,the field-theoretic approach of Martin, Siggia and 
RoseEZl (MSR). In this approach, the dynamic correlation 
and response functions 

C(r,t) = {u{v,t)u{Q,Q)) , 

[du][du] M(r,t)u(0,0)e-^1"'"1, (4.2a) 

R{v,t) = (7i(r,t)w(O,0)) , 

= y"[rfu][dM] u(r,t)u(O,0)e-^["'"l, (4.2b) 

are computed directly by integrating over the phonon and 
response fields u and {t, treated as independent stochas- 
tic fields with a statistical weight e~'^t"'"l imposed by 
the equation of motion, after integrating over the ther- 
mal noise C,{Y,t). The resulting effective "action" S is 
given hy S — Sq -\- Si, where 



+ iu{r,t)[-/dtu~ KV^u]^ (4.3) 

is the action of a pinning-free (harmonic) interface, and 
where Si = Sg + Sx, with 



/pFt 
drdt m(r, t) sin [pu{r, t) H ) (4.4) 



the contribution of the pinning potential and 
Sx[u,u\ ^ f drdt m(r,t) (Vu)^ 



(4.5) 



the contribution of the KPZ term to the nonlinearities in 
S. To study the renormalization of S[u, u] it is sufficient 
to work with the dynamic "partition function" 



[du][du] e-^["'^l. 



(4.6) 



required to remain fixed at unity under an RG coarse- 
graining procedure. The advantage of the MSR formal- 
ism is its close resemblance to the equilibrium statistical 
mechanics, that makes it a rather straightforward task 
to apply RG transformations and to derive recursion re- 
lations for the various parameters entering the equation 
of motion (4.1). Like the static momentum-shell RG, the 
dynamic RG procedure consists of three main steps : 

(i) Thinning of the degrees of freedom, whereby modes 
u{q), with q in an infinitesimal shell A/6 < q < A 
(6 — e'^^) are perturbatively (in Si) integrated out. 

(a) Rescaling of space variables according to r = 5r', 
so as to restore (for convenience) the ultraviolet cutoff to 
its original value A, and rescaling time variable according 
to t ^ t'b\ 

(in) Rescaling of fields, in order (for convenience) to 
keep the harmonic part of the action invariant under 
rescaling in (ii). 

We define "slow" {u<,u<} and "fast" fields {u>,u>} 



u(q,<) = u<{q,t) + u>(q,t) , 
M(q,<) = u<(q,i) +u>(q,t) , 



(4.7) 
(4.8) 



with momentum support in Fourier space in the inter- 
vals < q < A/b and A/b < q < A respectively, and 
perform a cumulant expansion of Z in terms of Si[u, u], 
considered as a perturbation. 



Z= I [du][du] e-So[n<,S<]^g-Si[«,f.]^^ 

[du][du] e-^"[""'""l-<^i)»>+^<^i)S>> , (4.9) 



where (• • •)o> denotes an average taken with the statis- 
tical weight So[u^ jU"^], and where the superscript c in 
('S'i)o> denotes a connected average. To first order in 
the pinning strength g, there is only one term in {Sg)^ , 
which renormalizes the dynamic action 5, namely 

{Sg)>^pgb-^r>V^^K^ 

X J drdt iu<{r,t)sm{pu<{r,t) + ^t), (4.10) 
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which physically arises from the suppression (from g 
to gb~^P q£ ^j^g effective pinning strength due to 

short-scale thermal fluctuations averaging away the peri- 
odic potential. In the above and throughout we will use 
= to indicates that only the leading term has been kept. 
Similarly, to first order in the KPZ coupling A, we have 
the following perturbative correction to the dynamic ac- 
tion S, 



{S. 



drdt iu (r, t) 



AttK 



(4.11) 



which quite clearly renormalizes the effective external 
force. 



Rescaling the space and time variables 



Y — hr' , 



as well as the conjugate field u(r, t) 
u<(r,t) = 6*u(r',t') , 



(4.12a) 
(4.12b) 



(4.13) 



while for convenience leaving M(r,t) unchanged in or- 
der to pjieserve the periodicity (27r/p) of the original 
problemlj we obtain the following lowest-order recursion 
relations : 



(7T)(6) 




(4.14a) 


lib) 


= 6'+* 7, 


(4.14b) 


K{b) 


= 6^+* K , 


(4.14c) 


gib) 




(4.14d) 


m 




(4.14e) 


{F/im 


^b^{Fh). 


(4.14f) 



The dynamic exponents z and x can be fixed by requiring 
that K and 7 be unchanged, to linear order in under 
the RG transformation. This leads to the following val- 
ues 

z = 2 , 

and to the following recursion relations for g and F 



^ = (2 



dF 



2F- 



4ttK 
ATA2 



)5, 



(4.15a) 
(4.15b) 



while the remaining quantities, K, 7, A, and tempera- 
ture T, remain unchanged and suffer no renormalization 
to first order in g and A. Similar considerations, with de- 
tails given in Appendix B, lead to the following recursion 
relations to second order in g and A : 



1 



SttK^ 167rA'3A4 1 + p 



(7T), (4.16a) 



m 

d9_ 

d£ 
dK 

m 

dX 

d£ 
dF 

m 



Tp'g' 1 - 

2\2 



IdnK^A^ (1 
Tp2 



1 , 



4:TTK 

Tp^g^ 2 
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3/2 - f 



lenK^A^ (1 -I 



J2)3 

K5) 



^ = 2F 



ATA2 Tp^g^ 



f 



AttK %'kK'^A^ 1 + /2 



(4.16b) 

(4.16c) 
(4.16d) 

(4.16e) 

(4.16f) 



wher e / = 
tion $A% 
the driven 
case (A = 
T{1) flows 



{j)F/2KA'^) is the dimensionless force of equa- 
. Note that, because of the lack of a FDT for 
system, in str ong co ntrast to the equilibrium 
F = 0) Eqs.( [4.16a| ) and ( [4.16b| ) imply that 
nontrivially according to 



dT _ . 



Tp^g^ 



P 



8nK^ SttK^A^ (l + PY 



T 



(4.17) 



Hence, T{£) is simply a measure of the strength of the 
white-noise component of the random force on the driven 
interface and is not associated with any equilibrium bath 
at a well-defined thermodynamic temperature. 

The recursion relations (4.16a)-( 4.16l| ) contain most 
(but not all, as discussed in Section V) of the informa- 
tion we need to investigate the properties of the system 



beyond the failing perturbative expansion of Section III. 
Before turning to their full analysis and to the study of 
the velocity-force characteristics, it is useful to see how 
the pr£j;ktusly derived static and equilibrium dynamic 
resultfoS are recovered. We do this in the following 
subsections. 



A. Analysis of the static limit 



The static model, Eq. (2.1), is characterized by two 



parameters K and g with the RG recursion relations re- 
ducing to the familiar Kosterlitz-Thouless form (clesived 
by these last authors in a dual. Coulomb gas formE^) 



d£ ^ 4ttK' 
dK _ Tp^g^ 
H ~ 87ri^2A4 ■ 



(4.18a) 
(4.18b) 



At small g, K{£) flows slowly, and the recursion rela- 
tion for g implies the existence a phase transition (called 
"roughening" in the context of crystal surfaces) at Tco = 
SttK/p^ (in the limit g ^ 0) between two phases distin- 
guished by the long scale {£ ^ 00) behavior of g{£). For 
T > Tco thermal fluctuations are strong enough to effec- 
tively average away the long- length scale effects of the pe- 
riodic pinning potential, which is therefore qualitatively 
unimportant for most (but not all) physical properties 
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of this so-called "rough" phase. At these high tempera- 
tures the surface is logarithmically rough and the effects 
of a weak periodic potential can be taken into account 
in a controlled perturbative expansion. In strong con- 
trast, for T < Tco, the effective strength of the periodic 
potential relative to that of the harmonic elastic energy 
grows on long length scales, leading to a breakdown of 
perturbation theory in g, no matter how weak its bare 
value might be. As a result, at long scales, the interface is 
pinned in this "smooth" phase, with bounded rms height 
fluctuations. 

It is instructive to recall some of the physics which 
follows from the above recursion relations. It is conve- 
nient to first rewrite the flow equations for dimensionless 
couplings g and rj 
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KA^ 
AttK ' 



(4.19a) 
(4.19b) 



that satisfy 



de 

drj 
dl 



1 2-2 

-^V 9 ■ 



(4.20) 
(4.21) 



These show that in equilibrium, the quantity rj which 
is the measure of the ratio of thermal (T) to elastic (K) 
energy, always flows to zero at long scales, indicating 
that the low-temperature smooth phase is controlled by 
a strong coupling zero-temperature fixed point. Near Tc, 
it is convenient to use a reduced temperature measured 
relative to the (noninteracting) Tco — SnK/p'^, 



= 7^-2, 

= 2(r/r,o - 1) 



with the flow equations simplifying to 
dg 

di^-^^' 
df _2 
11 = -' ■ 



(4.22a) 
(4.22b) 



(4.23a) 
(4.23b) 



Thes e can b e easily integrated by multiplying Eq. (4.23a) 
and (4.23b) by g and f, respectively, and taking the dif- 
ference of the two resulting equations. The result is that 
near T^o the flows are a family of hyperbolae 



1 . 5 




FIG. 4. Renormalization-group flow in the (f, g) plane. 
Temperature variation for an actual system occurs along the 
dashed line. On the high temperature side of the separatrix 
T = g (indicated as the thick line), the periodic pinning g 
renormalizes to zero and the interface is rough on long length 
scales. Below Tc (to the left of the critical separatrix) the RG 
flow run off to strong couphng g describing an interface that 
is smooth on long length scales. 



~2 -2 

9 ~T = c, 

labeled by a constant of integration 

2 



KA^ 



4^ 



- 2 



(4.24) 



(4.25) 



determined by the bare value of model parameters g and 
K. The resulting flows are illustrated in Fig. |4[ show- 
ing three distinct regions of behavior. In the high tem- 
perature region below the thick line (c < 0), pinning 
is irrelevant, and it therefore describes the rough phase, 
separated from the low-temperature smooth phase (the 
region above the thick line) by a critical line separatrix 
f = g. The latter therefore defines a true critical tem- 
perature given by 



1 



p^g 



V2A'A2 



(4.26) 



distinct from its g — > limit of T^o = SnK/p'^. Changing 
T corresponds to the variation of the dimensionless bare 
parameters along the dashed horizontal line indicated in 
Fig. I 

Above Tc, g{£) flows to zero and 



TR 



(4.27a) 
(4.27b) 



corresponding to the long-scale renormalized elastic con- 
stant 



Kii = Kie^oo), 

J-cD ^ 



(4.28a) 
(4.28b) 



It is comforting to find (using Eq. (4.25)) that AT/j re- 
duces to its bare value K at high temperatures. Using 
the fact that near, but above Tc 
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1 . 05- 
1 . 025- 
1 - 

Kr htl 

K 0.975- 

. 95- 

0. 925- 



1 2 3 4 5 

T Tco 

FIG. 5. Effective interface stiffness as a function of T /Tco 
for g — 0.1. In the smooth phase, Kn scales with the system 
size, and is effectively infinite. At T = , Kr takes the value 
Kr{T+) = K{1 +p^g/y/2KK^) with a universal ratio p^/Stt 
to the transition temperature Tc- Far above Tc, Kr goes to 
its bare value K. The dashed line indicates the location of 
Tc, which here is given by — 1.05Tco. 



~ -2f,T, 



(4.29a) 
(4.29b) 



with the true reduced temperature relative to the true 
(finite g) Tc given by 



2Tr\ T~Tr 



(4.30) 



and Tc = g = 2{Tc/Tco — 1), we find that in the limit 

Kn{T) ^K^(l~ ^g\T/Tc - . (4.31) 

This leads at Tc to a renormalized value of the elastic 
constant KjiiT^) that is enhanced relative to the bare 
value K and with the universal ratio to Tc given by 



Kr{T^c 



El 



(4.32) 



consistent with the analogoua_result first discovered in 
the context, of. the XF-modelO, related to our problem 
by dualityy^El. 

Below Tc, the relative pinning strength runs off to 
strong coupling and the interface is smooth on length 
scales longer than the correlation length that we calculate 
below. Because the RG flows are qualitatively very dif- 
ferent near and away from the two separatrices g = ±f , 
the value of this important length scale, that enters the 
velocity-force characteristics depends crucially on the dis- 
tance from Tc- In the critical region, defined by values 
of the bare parameters such that the weak-coupling (g) 
flow is near and roughly along either separatrix. 



l±~gf 



(4.33a) 
(4.33b) 



it is easy to show that the RG "time" 
coupling is given by 



to reach strong- 

(4.34a) 
(4.34b) 



Consequently the correlation length in this critical region 
is of familiar KT formO 



ae ' , 

„ga/|l-T/T.|V^ 



(4.35a) 
(4.35b) 



diverging extremely fast as T — s- , with a — \/2/g = 
{Tc/Tco — l)^"'^'''^ a nonuniversal constant. 

Outside this critical region, defined by r < —1, deep 
in the smooth phase, the flows are qualitatively different. 
At weak coupling gp^ ^ KA^ (the only regime where 
the perturbative RG analysis is valid) because f{£) grows 
weakly (additively) 



(4.36) 



grows exponentially fast, reaching strong-coupling at the 
low-T correlation length « ae^f given by 



\P 9 



(4.37a) 
(4.37b) 

(4.37c) 



On scales longer than the roughness correlation length 
the interface is smooth and is characterized by a strongly 
downward renormalized value of the pinning strength g ^ 
determined by the value of unrescaled coupling g{£ — 
log(^A)) at the scale of the correlation length. Near the 
transition 



t^t: 



-2a/\l-T/Ta\^^^ 



(4.38a) 
(4.38b) 



Deep in the smooth phase, for weak pinning, we instead 
flnd 



5fi ~ <7(A0"'' < 5, 

^ „2/(2~r,) 



T-^Tc 



(4.39a) 
(4.39b) 



which for weak g is also substantially reduced by thermal 
fluctuations. 

For strong pinning, fluctuations are unimportant and 
the correlation length reduces to the substantially shorter 
strong-coupling value — {K/ gp^Y^"^ determined by the 
bare model parameters. 



9 



B. Analysis of the equilibrium dynamics 

We now turn our attention to the equilibrium {F — 
A = 0) dynamics of the sine-Gordon interface, charac- 
terized by an additional model parameter, the friction 
coefficient 7, with the RG flow given by 



dl 1 ~2 



(4.40) 



Combining this with the recursion relation Eq. (4.18b), 
we find that the renormalized surface stiffness and 
friction coefhcient are related by 



7fl = 7 



(4.41) 



This together with the results of the previous subsection, 
show that the macroscopic linear mobility 7^^ is finitely 
renormalized in the rough phase, T > Tc and displays a 
square-root cusp approach to 7^^(T+) = "J^^iTco/TcY^^ 
as T ^ T+ 



1r 



1 



1/2 



(4.42) 



similar to the results of Petschek and ZippcliusS for 
the renormalized diffusion coefficient of the XY-model 
as T ^ T, 



KT- 



The effective friction coefficient ^{tj at scale e can be 
obtained by integrating the flow equation (4.40) 



7(£) =7 exp 



(4.43) 



Since below Tc, at weak coupling, g^{t)rj{i) grows with I, 
we find that the effective friction coefficient runs off to in- 
finity as ^ — > cx) suggesting a vanishing of the macroscopic 
linear mobility in the smooth phase. A more detailed 
analysis of the equilibrium weak-coupling flow equations 
for large ^ gives 



7(£) « 7 . 



exp 



exp 



|t|£ 



16(2-r,) 



T ^T: 



T<$:Tc 



(4.44) 



Such diverging friction coefficient can be physically inter- 
preted as activated creep dynamics over a pinning energy 
barrier that asymptotically grows with length scale, log- 
arithmically for T — > T~ and as a power-law for T <^ Tc- 
It is important to keep in mind tha t this growth of the 
friction coefficient ^{i) found in Eq. ( 4.44 ) extends only 
up to the stron g-coup ling length scale ^ = ae^' (£e for 



T T", Eq. ( tl.35b| ), and for T < T,, Eq. (|4.37c| )) 



since it was derived based on a renormalization group 
approach that is perturbative in g. In Section 0, we will 
look in more detail at the physics on scales longer than 
^, but we can already say at this point that (as we show 
in Section ^ even in this strong coupling regime the ef- 
fective 




FIG. 6. Effective linear mobility 7^ as a function of 
T/Tco in equilibrium {F — 0) for g — 0.1. Below the roughen- 
ing temperature at Tc, the mobility vanishes and the interface 
is pinned. 7^^(r) shows a square root cusp as T ^ , and 
goes to its bare value 7"^ for T ^ T^. The dashed line indi- 
cates the location of Tc, which here is given by Tc = l.OSTco. 

friction coefficient diverges. Consequently, we find that 
the interface linear (and in fact any order-n) mobility 
exhibits a nonuniversal jupip-discontinuity to zero across 
the roughening transitioiiE§c3, as illustrated in Fig. ^ 



V. NON-EQUILIBRIUM DYNAMICS AND THE 
VELOCITY-FORCE CHARACTERISTICS 



A. Weak-coupling regime 



We now turn to the full nonequilibrium problem, with 
the aim of deriving the velocity-force characteristics of an 
interface driven through a weak periodic potential, going 
beyond the failing (for T < T^) perturbative approach of 
section III, As long as the pinning remains weak, the long 



scale physics of the driven interface is contained in the 
renormalization group equations Eqs. 4.21-4.26, which 
when rewritten in terms of the dimensionless variables g, 
rj, f and the new KPZ coupling 



A 
pK 



(5.1) 



are given by 



c?7 

le 

dg_ 

d£ 
dl] 
dl 
d~X 

H 

dl 

di 



1 2^2-2 + 5/2 + 3/4 



(1 + /') 



2\3 



1 ~2/Cfl+5) 
S'^-^ (l+/2)3 ' 



(5.2) 
(5.3) 
(5.4) 

(5.5) 
(5.6) 
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The most striking effect of nonequilibrium dynamics is 
the breakdown of the FDT and as a result a nontrivial 
upward renormahzation (flow) of the effective "temper- 
ature" T(£) driven by the external force and the KPZ 
nonlinearity, reminiscent Q£.|naiiequilibrium "heating" in 
randomly-pinned systemsElfErQ. Consequently, even for 
T < Tc, for sufficiently strong drive the parameter 2—rj{£) 
determining the long-scale behavior of the periodic po- 
tential is driven negative, leading to the irrelevance of 
the pinning potential. Hence, as discussed in the Intro- 
duction, a finite external drive removes the qualitative 
distinction between the rough and smooth jphascs.and 
therefore rounds the roughening transitionJlaO'oH 

Here, we instead focus on the creep regime, where these 
particular nonequilibrium effects are unimportant. In 
this weak driving creep regime, we can ignore the KPZ 
nonlinearity and the most important role of F, as can 
be cl earl y seen even at the level of perturbation theory, 
Eq. (3.5), and from the equation of motion, is to intro- 
duce a new length scale ^ defined in Eq. ( p. 241 ). 
Beyond this nonequilibrium length scale the effects of the 
pinning potential and its ability to renormalize 7(f) and 
K{£) are suppressed, as it is averaged away on scales 
longer than (see for example the RG flow equations 
above and analysis below). Hence, for weak external 
drive F, the effective values of friction and interface stiff- 
ness parameters are given by "/{ip) and K(£f) renormal- 
ized by Gaussian equilibrium fluctuations up to length 
scale = e^^ • This therefore translates the strong £ de- 
pendence of "/{£) into strong F dependence of the macro- 
scopic mobility -f~^{F). Substit uting £_p, Eq. (3.24), in- 



side our equilibrium flow, Eqs. (4.44), and using 



v{F) ^ F/^{£p), 



(5.7) 



we immediately obtain the velocity-force characteristics, 
Eq. (1.2), quoted in the Introduction. 

applies as long as 



This prediction for v{F), Eq. (|lj), 
the relevant F probes length scales £_p on which the equi- 
librium weak-coupling flow equations remain valid. As 
discussed in the previous section, these flows in fact break 
down due to strong coupling effects (with g itself cut- 
ting off thermal Gau ssian flu ctuati ons) for length scales 
greater than ^, Eqs. ( 4. 351: ), ( 4.37c ). Hence, our predic- 
tions for v{F), Eq. ( |l.2D , remain valid only as long as 
^p < ^ (i.e., it is the external force and not the periodic 
potential itself that cuts off the Gaussian fluctuations), 
which translates into the condition F > F^,, with the 
crossover force given by equation (1.1) and in agree- 
ment with perturbation theory. 

To see this weak-coupling phenomenology emerge 
direc tly from our full nonequilibrium flow equations, 
Eqs.(5.2)-(p^), we integrate these equations, with A = 



and ignoring the nonequilibrium flow of T(£) (a valid ap- 
proximation in the F limit). We find for the renor- 
malized friction coefficient the following intermediate re- 
sult 



7fl(/) = 7 exp 



'd£r,{£)g\£) J_J^ 
{1 + P{£)y\ 



(5., 



Since at low drive and weak coupling, well-below Tc, 
r]{£), K{£), and T{£) grow slowly and f{£) and g{i) grow 
strongly according to 



m 



/e 



2£ 



(5.9a) 
(5.9b) 



it is quite clear from Eq. (5.8) that as long as the weak 
coupling flows remain valid, in the smooth phase the 
flows are automatically cut off when f{£) gets to be > 1 
leading to -^f discussed above. 

Su bsti tuting Eqs. (5_^) into the expression of jnif), 
Eq. ( |5.8| ), and integrating the resulting expression, we 
find 



7fl(/) = 7 



exp(ir;g2 A{^^J)) 



(5.10) 



with 



(2 - 77)/2 = (1 - T/Teo), and A(x, /) is the 



function given by 



2/4 V 2(2 



,Fi(2,2-x,3 



2fil-x)2F\2,l-x,2-x,--^) 



(5.11) 



gH{L 



g H{L 




FIG. 7. Behaviour of the pinning strength g (top) and 
of the friction coefflcient ^{£) (bottom) with length scale I 
for g = 0.1, T ~ O.STc and / ~ 1.7374 x 10"^ Here 
/, ~ 1.7373 X 10-^ 
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T TaO 

Characteristic force f*{T) as obtained from the 
numerical solution of the dynamic RG recursion relations 
(solid line) and from the perturbative estimate of equation 
( [3.22| ), (dashed line), for g = 0.1. The curve /.(T) delimits 
two very different physical regimes. Above this curve, the 
interface moves with uniform velocity. On the other hand, 
for / < /»(r), the interface moves through the nucleation of 
soliton excitations. 



where denotes a hypergeometric functionE3. When 
T < Tc {i.e. (—f) > ), ta king the limit of the function 
v4(— J, /) of equation ( 5.1C ) when / ^ leads to the fol- 



lowing expression for the long-scale inverse of nonlinear 
mobility 7/} 



7^(/)=7exp((^^./F)2(i-™), 



with 



(_nf_\ 

V 16(2-7;)/ 



(5.12) 



(5.13) 



in full agreement with earlier more qualitative discussion 
of the velocity-force characteristics in the intermediate 
regime of forces F_>_F* , and consistent with the per- 
turbative result (F22I) for 5 < 1. 

As F is lowered below F*, eventually the saturation 
of 7(£) breaks down and the flow behavior changes dra- 
matically as strong-coupling length scales (at which our 
weak-coupling RG solution is invalid) are probed. Study- 
ing the point at which this happens as a function of model 
parameters, allows us to extract the crossover value of , 
which we plot in Fig. ^. We flnd that there is a qualita- 
tive agree ment between the analytical prediction for 
Eq. (3.22), and our numerical analysis. 




(a) (b) 

FIG. 9. Schematic representation of the motion of a 
driven interface past the periodic pinning potential, (a) When 
T > Tc, or T < Tc and f > f*, the large fluctuations of the 
interface wash out the pinning potential on large length scales 
and the interface moves with a uniform velocity, (b) On the 
other hand, for T < Tc and f < f* the fluctuations of the 
interface are small; as a result, most of the interface is pinned 
at a given minimum of the pinning potential, and motion 
from one minimum to the next takes place through soliton 
excitations. 

the pinning barrier that scales like relative to the elas- 
tic energy. 

Instead, at low temperature the fluctuations are dom- 
inated by nontrivial saddle-point solutions (solitons) of 
H, Eq. ( |2.1[ ), with model parameters, Kn, gn, 7^ renor- 
malized by Gaussian fluctuations on weak-coupling scales 
L < ^. The dominant soliton excitation, illustrated in 
projection in Fig. ^ corresponds to a circular patch of 
radius i? > ^ of a nearly flat interface moving over to a 
neighboring minimum of the periodic potential, with an 
energy cost that clearly grows linearly with R 



Fsoiiton(i?) ~ pgniR 



(5.14) 



where gn (Eqs. (|3|), ( p9| )) and ^ (Eqs.( p35b| ) 



(4.37c)) strongly depend on the proximity to Tc. At zero- 
drive, the barrier to such solitonic motion simply diverges 
and linear mobility vanishes identically. A velocity-force 
characteristics in the weak drive F < {i.e., ^ < 
regime can be analyzed via scaling nucleation theoryEa. 
In this creep regime the interface is in near metastable 
equilibrium with F introducing a contribution 



B. Strong-coupling regime 



Ef{R) w -FnR^d 



(5.15) 



The weak-coupling behavior found in the p reviou s 
s ubsec tion only extends up to the scale ^, Eqs. (4.35b), 
(4.37c). Beyond this strong-coupling length, in the equi- 
librium model, the growth of g{£) and 7(^) is cutoff by the 
pinning potential, and an approach nonperturbative in g, 
where pinning is treated on equal footing with the elastic 
energy, is required. In this strong-coupling regime Gaus- 
sian interface fluctuations, considered so far, are strongly 
suppressed by 



to the effective free energy. Balancing Ep{R) against 
the soliton energy i?soiiton(-R) we find that solitons of size 
larger than a critical radius 



are unstable. In the F 
of solitons of size Rr ^ 



P9rC\ j_ 
27rd / F 



(5.16) 



-> limit, thermal activation rate 
1/F is quite clearly the hmiting 
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step for interface creep motion. We therefore find tha t 
the weak-couphng velocity-force characteristics, E g. (|l.2| ) 



crosses over, for F < F*, to that given by Eq. ( |1.3D m 
the Introduction, with 




A-nd 



(5.17) 



For vanishing temperature and strong bare pinning po- 
tential our asymptotic (for F < F^,) result for v{F) re- 
duces to that found in Refs. However, at large 
T < Tc and weak bare pinning we predict a strong 
thermal renormalization of the characteristic pinning en- 
ergy 



2 2 2 

p g a 



2 2 /-2 



(5.18) 



by thermal fluctuations on scales smaller than ^. 



An interesting and experimentally relevant general- 
ization of our results is a study of creep dynamics 
of a two-dimensional solid, driven through a one- or 
two-dimensional periodic potential, with applications to 
driven 2D colloidal crystals and vortices in supercon- 
ducting films. Despite of considerably different geom- 
etry, in equilibrium these systems display a pinned-to- 
floating solid transition closely related to the roughen- 
ing transition of 2D interfaces. However, new interest- 
ing ingredients arise. Some of the most imppr.tant ones 
are the nonequilibrium conventive-like termtOO'Ll, vec- 
tor phonon displacement and concomitant possible im- 
portance of dislocations. Combined with the consider- 
ably interesting behavior of the scalar sine-Gordon model 
studied here, we expect these to lead to even richer 
phenomenology. We expect that studies of these will 
shed considerable light on numerous experiments and 
simulationsH. 



VI. DISCUSSION AND CONCLUSIONS 

In this paper we have studied the creep dynamics of 
a two-dimensional interface driven through a periodic 
potential. Using dynamic renormalization group meth- 
ods and matching to strong coupling, we have calcu- 
lated the velocity-force characteristics across the inter- 
face roughening transition. Consistent with previous 
studies, we find a qualitative change across the transi- 
tion in the weak-drive velocity-force characteristics, with 
Ohmic transport for T > and a jump discontinuity 
in mobility across the transition. For T < Tc, in the 
asymptotic creep regime F <C F^{g,T) and for strong 
bare coupling, where transport is via soliton activation 
at all scales, we recover previously found results for the 
velocity-force characteristics v{F). However, for weak 
bare coupling and strong thermal fluctuations, we pre- 
dict an intermediate drive F > Ft,{g, T) nonlinear regime 
with a continuously varying (with T) exponent, which 
asymptotically crosses over to the strong coupling result 
with strongly thermally renormalized characteristic pin- 
ning barrier. Unfortunately, because the characteristic 
force F^:{g,T) that delineates between the intermediate 
drive regime and the strong-coupling regime coincides 
with the force marking the breakdown of the perturba- 
tive high-velocity expansion, we expect it to be difficult 
to observe this intermediate drive regime. 

The new physical picture which emerges from thje 
present study complements previously made predictionsu 
which were based on a more elementary perturbativp ja.p- 
proach, as well as known results for the mobilityOO at 
zero external drive. On the experimental side, the above 
picture may aiied some light on experiments such as those 
of Wolf et a/.cfl, who found that the growth velocity v of 
a surface of crystalline Helium 4 is strongly reduced at 
Tc from an Ohmic behavior v ^ F for T > Tc to an ex- 
tremely slow growth rate for T < Tc, a result which is 
usually explained in terms-pf an onset of creep motion 
via soliton like excitations^. 
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VII. APPENDIX : STATIC MOMENTUM SHELL 
RENORMALIZATION GROUP 



In this Appendix, we present technical details on the 
derivation of the renormalization group recursion rela- 
tions for the driven sine-Gordon model in 2 -I- 1 dimen- 



sions. For completeness, we shall begin in section VII A 
by showingiinssi the standard momentum shellEJ RG with 
hard cutoflOiEl can be applied to the static version of 
this problem before deriving the ful l dyna mic equations 
at nonzero external drive in section VII B . 



A. Static RG 



We decompose the field h{r) in the Hamiltonian (2.1) 
into high and low wavevector components 



such that 



hir) = h<{r) + h>ir) 
h<ir)^ r/i(q)e^'i-, 

Jq 



(7.1) 

(7.2) 
(7.3) 
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where /< = J^^' and /> = Xf/, denote in- 
tegration in momentum space over the ranges < |q| < 
A/6 and A/6 < |q| < A respectively. In terms of these 
high and low momentum fields, the equilibrium Hamilto- 
nian Ho[h] — J drif (V/i) can be written as the sum 

Ho[h] = Ho[h<] + Ho[h>] . 

We now want to integrate over the fast component 
/i^(r). To this end, we rewrite the partition function 
Z = J[dh] exp(-/3iJ) in the form (here /3 = l/T is the 
inverse temperature) 

/[d/i<][d/i>]e-''^«[''"l-'^^«["'l-'^^i[''"+''"l , 

^^,^<^^-PH„[h<] J ^^,^>^^-^H„[h>]^(}m[h<+h>] ^ 

\dh<]e->'Holh<]+PlnZ^^-0H,lh<+h>]^^^ , (7.4) 

where Zq = J[dh^] exp{—(3Ho[h^]), and where the sub- 
script (0 >) means that the average with respect to h-^ 



is performed with statistical weight exp(— l3 Ho\h ^]) / Zf^ . 
The term between angular brackets in Eq. ( [7.4[ ) is then 
approximated by a cumulant expansion 



^^-0Hi[h<+h> 



\ _ , _ {Hi)q> 1 2\c 

/o> ^ 2" +2r2^^'"> 



(7.5) 



where (^^i)§> denotes the second cumulant {{Hf — 
(i/i)^))p^. When re-exponentiated, equation ( ^^ ) leads 



to the result 

with the effective Hamiltonian 

Heff = {Hi)o> - — {H^)o> + ■ 



(7.6) 



(7.7) 



The averag£S_i|n Eq. (7.5) can be easily evaluated, with 
the resultsQEJ 



{Hi)o> = -.ge-^p'«'(°) J dr cos[p/i<(r)] , 
(i/2^;5^ = lg2g-p^G>(0) J rfrdr'[eP'^^("— "')-l]{cos(p(/i<(r) + ;i<(r')))+cos(p(/i<(r)-/i<(r')))} , (7.9) 



where (r — r' 



(/i>(r)/i^(r'))o> is the elastic prop- 



agator for fast fields (here Jq is the zeroth order Bessel 
function) 



G>(r 



T 



> g«q-(r-r') 



Tde 

2ttK 



Jo(A|r-r'|). 

(7.10) 



Given that G> (0) = Td£/27rK, we see that t he fir st order 
cumulant (7.8), after the rescalings ( [7.1^)-([7.19| ), leads 
straightforwardly to the recursion relation (4.18a) for the 
pinning strength g. On the other hand, since the "kernel" 



JC{r) 



pP'G>(r) 



1 



(7.11) 



takes appreciable values only for small values of its ar- 
gument, we see that the first term inside the integral 
in equation (|7.9| ) will contribute higher harmonic terms 
(~ cos (2ph{r)J) to the effective Hamiltonian, and hence 
we shall discard this term as irrelevant. In the second 
term, we shall make the approximation 

coa[p{h<{r)-h<{r'))] ~ 1 - ip^{h<{r) - h<{r')f , 

~ 1 - ^p\r - r'Ur-r')pdah<{r)dph<{r), (7.12) 

where, in going from the first to the second line, we made 
use of the Taylor expansion 

h<{r)-h<ir')^ir~r')^da.h<ir). 



Inserting ( 7.12| ) back into Eq. (7.9), we obtain the follow- 
ing approximation to the second cumulant (we here use 
the symbol = to indicate that we retain only the term 
correcting the stiffness K) 

X j dr'(r-r')^/C(r-r') . (7.13) 

Since G^ cx (i£ = In 6, we can expand the exponential in 
a Taylor series in G^ , 

gP^G>(r-r') _ ^ ^ ^2(.>(j. _ ^ i/(G(r - r'))' . 

(7.14) 

Now, the renormalization of K involves the integral 

j dr r2 /C(r) ^ J dr [e^'^^^"-) - l] . (7.15) 



Inserting the expansion ( 7.14 ) into this last expression, 
the first term gives a contribution 



J drr'G>{r)^-VlG>{q)\^^^ 



(7.16) 



which vanishes identically, since G^ (q) has support only 
on the shell A/b < q < A. The second term gives 
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J drr2/C(r) = dr r^{G>{r)y 

TV In 6 



7ri^2A4 



(7.17) 



Thus, we obtain for the second cumulant ( 7.13 ) the fol- 
lowing expression 



{Hi 



o> 



I dr^{Vh<ir)y 



2T 8t:K^A^ 
We now perform the following rescalings 



momentum regions, respectively. Using the fact that 



u{q,uj) 



■(q, w) + u>(q, w), it is not difficult to verify 



that the free part of the action decomposes into two 
diagonal pieces Sq and Sq depending only on M^(q, a;) 
and u> (q, lu) respectively 



So[u,u] = So\u' 



(7.24) 



As we did in the static RG, in order to be able to in- 
tegrate out the fast component of the field u{r,t), we 
rewrite the generating functional Z in the form 



e'r' 



h{r') 



(7.18) 
(7.19) 



so as to restore the ultraviolet cut-off back to A. Because 
the pinning potential is a periodic function, it is conve- 
nient (although not necessary) to set the arbitrary field 
dimension x to zero, thereby preserving the period 2tt/p 
of the original problem under RG transformations. Un- 
der such a transformation, the resulting effective Hamil- 
tonian can be cast into its original form with effective 
^-dependent parameters K{£) and g{£) such that 



[du<][du^ 



]+lnZ> I -Silu<+u>, 



>0> 



where Zq 



exp(— S'o[u^, u^]), and where 
denotes statistical averaging with statistical 
weig ht e-^«["^'"^l. The perturbative correction to the 
dynamic action can therefore be expressed in terms of a 
cumulant expansion 



= 1 - 



1 



i;o> 



(SI) 



o> 



(7.25) 



K{t) = K - 
or, in differential form 



d£ 



^ = (2- 
di 



Tp^ 



dK 

m 



(7.20a) 
(7.20b) 

(7.21a) 
(7.21b) 



B. Dynamic RG 

We now turn our attentio n to t he deriv ation of the dy- 
namic RG flow equations ( 4.16a )-( 4.16i ) for the driven 
sine-Gordon model. As we did in the static case, we de- 
fine the following low and high momentum components 
h<{r,t) and h>{r,t) : 



'(r,i) 



> 



/i(q,w)e'('l--'"*) 
/i(q,w) e'(i-'-'^*) 



(7.22) 
(7.23) 



where, here and in what follows, f on integrals stands 

for J ^^J^a I and the superscripts > and < indicate inte- 
gration over the high {A/b < q < A) and low (0 < g < A) 



ASc 



\u<u] 



Reexponentiation of this expansion allows us to define 
the effective action 



0> 



(7.26) 



from which we can derive dynamic RG flows for the pa- 
rameters of the original equation of motion. This proce- 
dure, to first order in the pinning strength has already 
been shown in the text. Here we are therefore only going 



to consider the second ordei; 
tion S. In fact, it turns outi 



sorrection to the original ac- 
3 that the only perturbative 



corrections to S to second order in perturbation theory 
come from the cumulants —\{S'^)oy and — ^('S'^)o>) *-e. 
we need not consider the cross term —{SgS\)oy which 
does not provide any perturbative corrections to the ac- 
tion. In the following, we shall only show how we com- 
pute the perturbative corrections arising from the sine- 
Gordon perturbation — i(S'g)o>, the unique term arising 
from-i(52)o> 



ASx{lT) = 



dvdtu\ (r, t) 



TX^ di-] 



(7.27) 



having been repeatedly derived in the litera turdO 'ElB. 
Taking the Gaussian averages in equation (7.25) leads 
to the following expression of the second cumulant 



ASg[u, u 



dvdt dv'dt' U<{Y,t)u< {y' ,t')t{Y - t' ,t 



t')cos [p{u<{y, t) 



nr',t')) 



dYdt 



dY'dt' iu<{Y, t)/C(r - r', t - t') sin [p(u<(r, t) - u<(r', t')) + —{t- t')] 



(7.28) 
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by 



Here the dynamic kernels JC{r,t) and JC{r,t) are given CQ{r,t) = ([w^(r,t) — u-'(0, 0)]^) are the response and 

correlation functions, respectively, and where the corre- 
lator G^{r,t) = (u>(r, i)u>(0, 0))o> is given by 



^(r,i) = -[l-cosh(p2G>(r,i))] 

-sinh(p2Go>(r,t)) , (7.29) 
JC{r,t)=e-ip'^oi-'t)R>(^r,t) , (7.30) 



G>{r,t)^2^T 



cos[q ■ r — ujt] 



(7.31) 



We now decompose the sine and cosine in the integrand 



where Ro{r,t) = J^^^e '(i"" "*V(«7^ + ^9^) and on the rhs of equation ( [7.2^ ) according to 



cos [p{u< - it'<) + ^{t - t')] = sin [p(u< - u'^)] cos [^{t - t')] + cos [p{u< - u'<)] sin [^{t - t')] , (7.32a) 



sin [p(u< - O + ^(t - t')] = sin [p(u< - it'<)] cos [^{t - i')] + cos [p(w< - u<)] sin [^{t - t')] • (7.32b) 

The kernels /C(r — r', t — t') and /C(r — r', i — f') being short ranged both in space and time, we see that the major 
contribution to the action (7.28) comes from the regions r ~ r' and t ~ t' where [it^(r, t) — u<(r',t')] is small. We 
therefore shall approximate 



sin[piu<ir,t)-u<ir\t'))] ^ p {u< ir,t) - u< (r^t')) , 

cos [p(u<(r, t) - u<(r', t'))] ^l-\p' («<(r, t) - u<{v', t')f , 



and 



u<{v,t) - u<{v', t') = {t- t') dtu< + (r - r')a d^u< + ^{r - r')a(r - r')pdo.dpu< 
upon which we obtain the following expression for the second cumulant —^{Sf)Qy : 



(7.33) 
(7.34) 

(7.35) 



where 



AS[u,u] = --{SYoy[u,u] = A5(7T) + ASi^f) + ASiK) + A5(A) + ASiF) 



fPF , 



AS{-fT) ^ ASxijT) + -p^g^ J drdt u<{Y,t)u<{r,t) J dr'dt' K-ir - r' ,t - t') cos[ — {t - t')] , 
AS'(7) = ipVy" drdtiu<{r,t)[dtu<ir,t)] J dr'dt' {t - t') IC{r - r' ,t ^ t') cos[^{t - t')] , 
AS{K) = ^p^g^ J drdtiu<{r,t)[-V^u<{r,t)] J dv' dt' {r ^ v'f JC[v - v' ,t ^ t') cos[^{t - t')] , 
AS{\) ^^p^'g^ j dvdt iu<{v,t)[ - {Vu<[v,t)f] J dv'dt' (r-r')'/C(r-r',i-t')sin[^(t-t')] 
AS{F) = ^p^g^ J drdtiu<{r,t) J dr'dt' JC{r - r' ,t - t')sin[—{t - t')] . 



(7.36) 

(7.37a) 
(7.37b) 
(7.37c) 
(7.37d) 
(7.37e) 



Here we pause a moment to indicate that if we use the 
complete expression of the kernel /C(r, t) 



/C(r,i) =e~5P'Co>(r,t) R>(r,t) 



(7.38) 



into equation ( 7.37e ) and let b —> oo, then we obtain 
from equation ( 7.376 ) above the following expression for 
the friction force Ffr due to the pinning potential to or- 
der g^, 

Ffr ^Ip^9^ J rfr'dt'e-if'^°('-'-''*-*'' x 



X i?o(r-r',f-i')sin[— (i-t')] , (7-39) 
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which leads directly to the perturbative result (3^) of 
the text. 

We now g o back to our dynamic RG recursion rela- 
tions, ( 7.37a )-( 7.37e ). In the dynamic kernels of equa- 
tions ( 7.29| )-( 7.3Cl| ), we expand 



^(r,t) = -p'G>{r,t) - ^p\G>{r,t)Y , (7.40a) 
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/C(r, t) = R> (r, i) - - p^C> (r, i)i?> (r, t) , (7.40b) 

and keep only the second term on the rhs of the above 
equations (the first term gives a vanishing contribution, 
for reason s wh ich are identical to those explained after 
equatio n (|7.16| ) of appendix A). Now, from equations 
( 7.37b| )-( [7.37eD , we see that the perturbative corrections 
to the bare parameters of the theory are given by the 
flows 



dilT) 




TA2 


M 


pert 








1 4 2 


m 


pert 




dK 




1 4 2 


~M 


pert 


IP 9 


dX 






H 


pert 




dF 






'dl 


pert 





2 2 

P 9 



pF 

dr dt /C(r, t) cos ( — t) , 

7 



pF 

dr dt tfCiv, t) cos ( — t) , 

7 

vF 

drdtr^lCir,t) cos(—t) 

7 

j drdtr^JC{r,t) sin (— t) 
pF 

dr dt /C(r, t) sin ( — i) . 



Using equations (7.40a)-(7.401: ), the above recursion re- 
lations become 



dj 

m 

dg_ 
d£ 
dK 
~d£ 
dX 

n 

dF 



Tp'g' 



1 



UttKS ' WttK^A^ 1 + /2 
Tp^g^ 1 - f 



IGttK^A'^ (1 + /2)2 
. Tp^ 



7 



9 



AttK' 
Tp^g^ 2-3/2-/4 



167rif 2A4 (1 + /2)3 
W /(/^ + 5) 

ArA2 rp5c,2 J 



47rii: 8t:K^A^ 1 + /2 



(7.41a) 
(7.41b) 

(7.41c) 
(7.41d) 

(7.41e) 

(7.41f) 



On the other hand, we know from equations (4.12a)- 
(4.13) that the rescaling of fields and space and time 
variables produces the recursion relations 



d£ 
dF 

d£ 



dj I 

di ''■'^^'^ 



dK. _d\. _ ^ 



^ =2F 



(7.41g) 



Using the recursion relations above along with the fact 
that, in a renormalization group transformation. 



d _ d 
dl~d£ 



+ — 

pert d£ 



(7.42) 
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